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The "complementary" Ansatz, Tr(Mi,)=0, where Mu is the prediagonal neutrino 
p^ ' mass matrix, seems a plausible approximation for capturing in a self contained way 

\-^ ' some of the content of Grand Unification. We study its consequences in the form 

of relations between the neutrino masses and CP violation phases. 



1. Introduction 



A favorite topic of discussion at the MRST meetings over many years has 
been the prediction of quark and lepton masses and mixings. Usually, an 
attempt is made to predict everything at once based on a suitable guess 
(Ansatz). Here we discuss a "complementary" Ansatz which adds just a 
little information to the system. This is ^'^''^ , for the symmetric pre- 
diagonal neutrino mass matrix, M^: 

Tr{M,) = (1) 
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If, at first, leptonic CP violation is neglected so that M^ is diagonalized by 
a real orthogonal matrix, Eq. (1) yields simply 

mi + 1712 + m^ — 0, (2) 

where the neutrino masses, rrii can be taken here to be either positive or 
negative. Now a very recent analysis * of solar, atmospheric, reactor and 
accelerator neutrino oscillation data (but neglecting LSND) gives the best 
fit: 

ml -m\ = 6.9 x lO^^e^^ 
Wi -ml\ = 2.6x IQ-^eV'^. (3) 

Together, Eqs. (2) and (3) provide three equations for three unknowns. 
There are two essentially different types of solutions. Type I is characterized 
by ItosI being largest: 

mi = 0.0291 eV, ma = 0.0302 eV, mg == -0.0593 eV, (4) 

while type II has |?n3| smallest: 

mi = 0.0503 eV, mz = -0.0510 eV, mg = 0.00068 eV. (5) 

Here we will, following ref.^ (see this for further references), discuss the 
situation when CP violation effects are included and give an application to 
leptogenesis. 

2. Plausibility argument for Ansatz 

SO(IO) grand unified theories have the nice feature that they contain a 
complete fermion generation in a single irreducible representation of the 
group. There are three Higgs irreducible representations which can directly 
contribute to tree level fermion masses via the Yukawa sector: the 10, the 
120 and the 126. In principle any number of each is allowed. For every 
Higgs field there is a 3 x3 matrix of unknown coupling constants. The 
fermion mass matrices are linear combinations of these matrices. Clearly a 
very large number of different models can be envisioned. 
We start from the "kincmatical" relation: 

rr(M_i/3 - rAf_i) ex Tr{Mo,LiGHT), (6) 

which holds when any number of lO's and 120's are present but only a 
single 126. Here the subscript on the mass matrix indicates the electric 
charge of the fermion. The quantity r « 3 takes account of running masses 
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from the GUT scale to about 1 GcV. Under the same conditions one also 
has: 

Mq.heavy oc Mo^light- (7) 

Note that the physical light neutrino mass matrix, Mi, is given by the well 
known formula: 

Ml, « Mq, LIGHT - Mdirac^o,heavy^dirac- (8) 

The initial assumption we shall make is that the second, "see-saw" term in 
Eq. (8) is small compared to Mq^jq^^t- 

Now, it has been known for a long time that the quark mixing matrix is 
of the form diag{l, 1, 1) + 0(e). Thus it was very surprising when analysis 
of neutrino oscillation observations showed that the lepton mixing matrix is 
not at all close to the unit matrix but rather has large (12) and (23) mixing 
elements. In a GUT framework, this suggests that a first approximation to 
the prediagonal mass matrices might be to take the charged fermion mass 
matrices to be diagonal while the neutrino mass matrix would presumably 
differ drastically from the diagonal form. As examples we would set M_i « 
diag{me,m^,mr) and M_i/3 « diaglrndTms^mi,). Substituting these into 
the left hand side of Eq. (6) shows it to be about {mi, — Smi-), which is 
about zero. Hence the right hand side should also be about zero as should 
Tr{M„) in the non-seesaw dominance case. Although this is clearly an 
approximation, it seems likely to be close to the physical situation in the 
same sense as nib ~ SttIt-. Of course, the approximation gets better as the 
mixing matrices needed to bi-diagonalizc the charged lepton mass matrix 
get closer to the unit matrix. For simplicity, in what follows we shall also 
approximate these matrices to equal the unit matrix. 

3. Parameterized Ansatz equation 

The prediagonal neutrino mass matrix may be brought to diagonal form by 
a transformation: 

U M^U = M^ — diag{mi,in2,m3), (9) 

where C/ is a unitary matrix. M^, is a symmetric but complex matrix 
which has in general 12 real parameters. This equals the sum of the three 
parameters from rrii and the nine parameters from U. Now the observable 
lepton mixing matrix, K is given ^ hy K = fl'^U, where we have just agreed 
to approximate the charged lepton diagonalizing matrix factor, ft to be 
essentially the unit matrix (or alternatively we could choose to work in 
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a basis where the charged leptons are diagonal). Thus we replace U by 
the observable matrix K. K is parameterized in a conventional way as 

K = Kexp^Q^^iT), where ujq{t) ~ diag{e"^ , e^'^'^ , e'^'^^) with Ti + T2 + T3 = 
and: 



K 



C12C13 S12C13 siae 

-S12C23 - Ci2Si3S23e*'^ C12C23 - Si2Sl3S23e*'' C13S23 
S12S23 - Ci2Sl3C23e*'^ -C12S23 - Si2Si3C23e*'^ C13C23 



(10) 



where s^ — sinOij and Cij — cosOij. As written, the matrix K is parame- 
terized by three angles and three independent CP violating phases. To get 
the most general K three more phases are required; these can be inserted 
for example ^ by multiplying K on the left by ojo {(t) . However these phases 
can always be canceled by rephasing the (diagonal) charged lepton fields 
which sit to the left of K. Thus even if the phases, Ci were included in 
U there would always be an allowed choice of charged lepton phases which 
would cancel their effect when we get restrictions (as we shall) on the phys- 
ical K. Taking this into account the Ansatz reads in terms of physical 
quantities: 

Tr(M,X-ipX:,p^o(2T)) = 0. (11) 

In more detail it is: 

mie^'^i [1 - 2i(ci2Si3) 



TO2e ^"^^ [1 — 2i(si2Si3) sinSe 



m3e-"'- [1 + 2i{si3)^sinSe'^] = 0. (12) 

In this equation we can choose the diagonal masses mi, 1112,1713 to be real 
positive. The mixing angles are known from the best fit ^, 



'12 



= 0.30, s^3 == 0.50, Si3 = 0.003, (13) 



wherein the first two have about 25 per cent uncertainty while the third is 
just known to be small. 

Together Eqs. (3) and (12) are now seen to provide 4 real equations for 
the 6 unknowns: mi, m,2,m3, S, ti,T2- Thus further assumptions are needed 
to make some predictions. We already saw that assuming all the CP phases 
to vanish gives three equations for three unknowns. If we assume only the 
"conventional" CP phase 5 not to vanish there are four equations for four 
unknowns, with the results described in ^ . This case would become trivial 
in the limit where 5^3 is assumed to vanish. There is no reason to expect 
it to vanish exactly but, considering our lack of knowledge, that seems to 
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be also an interesting assumption to investigate. According to Eq. (12) it 
yields the same result as setting (5 = 0. Then we have 4 equations for 5 
unknowns and can get a one parameter family of solutions. 

4. Family of Majorana phases 

In this case, Eq. (12) takes the form 

mie^'^i + mae^*^^ + m^e^"-^ = 0, (14) 

which corresponds, as illustrated in Fig. 1, to a triangle in the complex 
plane. 




Figure 1. Vector triangle representing Eq. (14). 

We proceed to obtain the family of solutions by assuming a value for 7713, 
getting TTii and m2 from Eq. (3) and finally by solving for the two interior 
angles /xi and ^2 using trigonometry. Interesting CP violation quantities 
turn out to be: 

sin[2{Ti - T2)] = -sin{^i + ^2), 

sin[2{Ti — Ts)] = sinfi2, 

sin[2{T2 - T3)] = -sinfii, (15) 

The criterion for the existence of CP violation effects is the area of the 
triangle being different from zero. We may express the area as: 

Area = - {[{mi + 7712)^ — m3][TO3 — (mi — 1712)'^]) . (16) 
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Now we see that the vanishing of the first factor corresponds to the type 
I real solution while the vanishing of the second factor corresponds to the 
type II real solution. Furthermore, for a solution to exist, the argument of 
the square root should be positive. With the second factor, that establishes 
the minimum allowed value of TO3 while with the first factor, that establishes 
the minimum value of TO3 which allows a type I solution. 

The table below shows a "panorama" of solutions decreasing from m^ = 
0.3 eV, (which is about the highest value compatible with the cosmology 
bound ® that the sum of the neutrino masses be less than about 1 eV) to 
the lowest value imposed by the model. In the type I solutions m^ is the 
largest mass while in the type II solutions 7713 is the smallest mass. For 
each value of TO3, the values of the model predictions for mi and TO2 as well 
as the internal angles /xi and /X2 are given. The model prediction for the 
neutrinoless double beta decay quantity |?Tiee| is next shown. Finally, the 
last column shows the estimated lepton asymmetries due to the decays of 
the heavy neutrinos. Note that the reversed sign of lepton asymmetry is 
also possible. 





mi, 7712,1713 in eV 


/ii,^2 rad. 


\mee\ 


ei, £2, £3 


I 


.2955, .2956, .300 


1.038, 1.039 


.185 


.342, .433, .017 


II 


.3042, .3043, .300 


1.055,1.056 


.187 


.330, .426, -.0172 


I 


.0856, .0860, .100 


0.946, 0.952 


.058 


.138, .060, .00137 


II 


.1119, .1122, .100 


1.106, 1.111 


.065 


.194, .088,-. 0024 


I 


.0305, .0316, .060 


0.258, 0.268 


.030 


.00982, .00422, .00004 


II 


.0783, .0787, .060 


1.172, 1.187 


.043 


.094, .041,-. 0011 


I 


.0291, .0302, 


.000552, 


.030 


1.96 xl0-^ 




0.0592715649 


0.000574 




.84 xl0-^ .71 xlO-^ 


II 


.0774, .0782, 
0.0592715649 


1.174, 1.188 


.042 


.047, .020, -.0011 


II 


.0643, .0648, .040 


1.243, 1.268 


.033 


.052, .023,-. 000681 


II 


.0541, .0548, .020 


1.355. 1.442 


.024 


.018, .0078,-. 000335 


II 


.0506, .0512, .005 


1.386, 1.658 


.021 


.0057, .0025,-.0000824 


II 


.0503, .0510, .001 


0.814, 2.313 


.021 


.00073, .00031, 
-.0000122 


II 


0.0503, 0.0510, 


0.051361, 


.021 


.0000348, .0000150, 




0.0006996 


3.089536 




-0.601 xlO-6 
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Notice that when m^ is decreased to about 0.0593 eV, we get to the 
real type I case (no CP violation) . Below this value of ms only the type II 
solutions exist. At TO3 about 7 x 10^^ eV, we get the real type II case and 
no solutions exist for 1713 below this value. In the TO3 regions just above the 
two real cases we can evidently tune the CP violation phases continuously 
to be as small as desired. 

An interesting application of the model is to neutrinoless double beta 
decay (for example, the decay ^^Ge -^ ^^5e + e^ + e^). The current 
experimental bound on the amplitude factor ^ is: \mee\ < (0.35 -^ 1.30) 
eV, where 

Imeel = J2 \MKe.pu)^e~^'^^\. (17) 

From the table we see that the predicted values of \mee\ are typically about 
one order of magnitude below the experimental bound. Furthermore, the 
predicted values do not vary drastically for 7713 less than about 0.1 eV. 

5. Estimate for leptogenesis 

An intriguing possibility for learning more about leptonic CP violation is 
the study of the proposed leptogenesis mechanism for generation of the 
present baryon number asymmetry of the universe. According to this 
scheme, the lepton number violating decays of the heavy neutrinos at a 
high temperature (early universe) establish a lepton asymmetry which gets 
converted as the universe cools, through a (B+L) violating but (B-L) con- 
serving "sphaleron" interaction to a baryon asymmetry. References and a 
rough estimate in the present framework are given in ref. ^. According to 
Eq. (7) the heavy neutrino masses are supposed to be proportional to the 
light ones here and have the same diagonalizing matrix, U. The effective 
term for calculating the heavy neutrino decays at very high temperature is 

Lyukawa = -Y^ Uh^^^Nj + H.c, (18) 

where 

'^^1- <</,o>r' • ^ ' 

The quantities needed for the calculation are the matrix products {h^h)ij; 
it is thus seen that the effect of a diagonal matrix of phases multiplying 
Kexpij on the left would cancel out. The lepton CP asymmetry e^, due to 
the decay of the ith heavy neutrino, is defined as the ratio of decay widths: 

T{N^ ^ L + $) - r(iv, ^ L + i>) 



T(N, ^ L + $) + r(7V, ^ L + $) 



(20) 
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In this equation L + $ stands for all lepton- Higgs pairs of the types e 



T 

and Vj+(t)^ . This is an effect which violates C and CP conservation, in agree- 
ment with the requirement of Sakharov. The numerical values of the e^, 
which depend on the ratios of heavy neutrino masses rather than their ab- 
solute values, are displayed in the last column of the table. Notice that Eq. 
(18) represents the same term which generates the Dirac mass, Mdiuac 
in Eq. (8). Since our motivation for the Ansatz assumes dominance of the 
non seesaw term, this feature requires ^ the heavy neutrino mass scale to 
be suitably large. This scale plays a role in the estimation of the present 
baryon to photon ratio, tjb of the universe which is obtained by convoluting 
the ei with factors obtained by solving the Boltzmann evolution equations 
for the (B-L) asymmetry. It turns out ^ that for typical values of the pa- 
rameter, 7713 in the table, tjb is considerably larger than its experimental 
value ® of about 6.5 x 10^^°. Thus agreement with experiment requires 
tuning close to the two real type solutions; the correct order of magnitude 
is obtained when either m^ « 0.059 eV (type I) or nis « 0.005 eV (type 
II). 

We thank D. Black, A. H. Fariborz, C. Macesanu , M. Trodden and D. 
Schechter for their help. The work of S.N is supported by National Science 
foundation grant No. PIIY-0099544. The work of J.S. is supported in part 
by the U. S. DOE under Contract no. DE-FG-02-85ER 40231. 

Finally, we are pleased to take this opportunity to express our wishes 
for Happy Birthday, Good Health and Continued Important Contributions 
to Physics to PAT O'DONNELL and to HARRY LAM. 
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